Soil resistances for laterally loaded rigid piles in multi-layered elastic soil

ABSTRACT

Short stubby piles like monopiles and large diameter drilled shafts undergo rigid body translation and rotation when subjected to a lateral force and/or a moment at the head. A method of analysis for these piles embedded in multi-layered elastic soil is developed using the variational principles of mechanics. Using this analysis, the soil resistance against pile movement can be rigorously related to the soil elastic constants, and the pile head displacement and rotation can be quickly calculated. The equilibrium equations for pile and soil displacements are obtained using the principle of virtual work and solved using an iterative algorithm. Pile responses obtained from the analysis match well with those obtained from three-dimensional finite element analyses in which the same inputs of loads, geometry, and material properties are given. Based on the new analysis, fitted equations for soil resistance parameters are developed, which can be used to directly calculate the pile head displacement and rotation without the use of the iterative algorithm. Numerical examples are provided that demonstrate how the method can be used to analyse practical problems.     

基于ROMS和4DVAR的沿轨与网格化SSH数据同化效果评价

Remote sensing products are significant in the data assimilation of an ocean model. Considering the resolution and space coverage of different remote sensing data, two types of sea surface height(SSH) product are employed in the assimilation, including the gridded products from AVISO and the original along-track observations used in the generation. To explore their impact on the assimilation results, an experiment focus on the South China Sea(SCS) is conducted based on the Regional Ocean Modeling System(ROMS) and the four-dimensional variational data assimilation(4 DVAR) technology. The comparison with EN4 data set and Argo profile indicates that, the along-track SSH assimilation result presents to be more accurate than the gridded SSH assimilation, because some noises may have been introduced in the merging process. Moreover, the mesoscale eddy detection capability of the assimilation results is analyzed by a vector geometry–based algorithm. It is verified that, the assimilation of the gridded SSH shows superiority in describing the eddy's characteristics, since the complete structure of the ocean surface has been reconstructed by the original data merging.     

Groundwater contaminant source identification based on iterative local update ensemble smoother

Identification of the location and intensity of groundwater pollution source contributes to the effect of pollution remediation, and is called groundwater contaminant source identifcation. This is a kind of typical groundwater inverse problem, and the solution is usually ill-posed. Especially considering the spatial variability of hydraulic conductivity field, the identification process is more challenging. In this paper, the solution framework of groundwater contaminant source identification is composed with groundwater pollutant transport model (MT3DMS) and a data assimilation method (Iterative local update ensemble smoother, ILUES). In addition, Karhunen-Loève expansion technique is adopted as a PCA method to realize dimension reduction. In practical problems, the geostatistical method is usually used to characterize the hydraulic conductivity feld, and only the contaminant source information is inversely calculated in the identifcation process. In this study, the identification of contaminant source information under Kriging K-field is compared with simultaneous identification of source information and K-field. The results indicate that it is necessary to carry out simultaneous identification under heterogeneous site, and ILUES has good performance in solving high-dimensional parameter inversion problems.     

非饱和土中水平入渗方程显式求解

在求解非饱和态土中水分入渗问题时,水力函数是体积含水率或者吸力的函数,致使其控制方程呈现出强非线性的特征,进而使得其求解变得十分困难。基于水分在土体介质中流动耗时取极值路径的选择这一假定,引入时间泛函,基于变分法原理将水平入渗问题转化为泛函极值问题。通过求解Euler–Lagrange方程,结合边界条件,得到非线性瞬态水平入渗问题的显式解析解。结合Brooks-Corey型水力函数,显式地求解出该类型非饱和态土的体积含水率发展分布规律。通过计算4种不同类型土体的水平入渗规律,将求解结果与已有结果以及数值结果进行对比,验证了该方法的有效性。结果表明:体积含水率分布与位置距离和湿润峰距离比值呈幂函数关系,指数取决于土-水特征曲线的形状参数;初始条件与边界条件会对体积含水率分布造成不同程度的影响。     

基于WACCM+DART的临近空间SABER和MLS臭氧观测同化试验研究

本研究在WACCM+DART（Whole Atmosphere Community Climate Model,Data Assimilation Research Test-Bed）临近空间资料同化预报系统中加入SABER（Sounding of the Atmosphere using Broadband Emission Radiometry）和MLS（Microwave Limb Sounder）臭氧观测同化接口,并以2016年2月一次平流层爆发性增温（SSW）过程为模拟个例进行了SABER和MLS臭氧观测同化试验,得出以下结论:同化SABER和MLS臭氧体积浓度观测得出的WACCM+DART臭氧分析场能够较真实反映SSW期间北极上空平流层臭氧廓线随时间的演变特征,且与ERA5（Fifth Generation of ECMWF Reanalyses）再分析资料描述的臭氧变化特征具有很好的一致性;基于SABER和MLS臭氧观测的WACCM臭氧6 h预报检验表明同化臭氧观测对臭氧分析和预报误差的改善效果主要体现在南半球高纬平流层和北半球中高纬平流层中上层-中间层底部;基于ERA5再分析资料的WACCM+DART分析场检验表明同化SABER和MLS臭氧体积浓度资料可在提高北半球高纬地区上平流层-中间层底部臭氧场分析质量的同时减小该地区上平流层-中间层底部温度场和中间层底部纬向风场的分析误差;基于MLS臭氧资料的臭氧中期预报检验表明相对控制试验同化SABER和MLS臭氧体积浓度资料能更好改善0~5 d下平流层和中间层底部臭氧的预报效果。     

适应数值模式垂直坐标的GPS掩星反演资料稀疏化方案

GPS掩星反演大气温湿资料具有高垂直分辨率、高精度、受云和降水影响小等优点,针对GRAPES同化预报系统,发展设计了一种既考虑预报模式高度-地形追随垂直坐标不均匀分层特点,又结合掩星反演资料特性的新适应性垂直稀疏化方案。通过个例试验和批量试验,探索了该适应性稀疏化方案对分析预报质量的影响。试验结果表明:选取合适的稀疏参数,新稀疏化方案的GPS掩星反演资料对背景场的调整更加有效,分析场质量更高;位势高度、比湿、温度和风场等预报场的均方根误差均更小;适应性稀疏化方案对改善台风路径预报具有积极作用。批量试验则进一步证实了适应性稀疏化方案对分析场质量有明显的改善作用。     

高时空分辨率三维风场在强对流天气临近预报中的融合应用研究

以业务应用为目标,开展高时、空分辨率三维风场在强对流天气临近预报中的融合应用研究。运用北京快速更新多尺度分析和预报系统集成子系统（RMAPS-IN,Rapid-refresh Multi-scale Analysis and Prediction System-Integration）,对雷达四维变分分析系统（VDRAS）30 min临近预报的高时、空分辨率三维风场作为数据源与自动气象站风场观测进行快速融合处理。结果表明,以VDRAS临近预报风场取代数值模式预报场作为融合初猜场后形成的分析结果对于风场有明显的改善:（1）长时间序列客观检验结果表明,地面10 m高风场U/V分量绝对误差分别为0.05和0.06 m/s。实时融合对未来预报的影响随着预报时效的延长,U/V分量的绝对误差不断增大。（2）对于11个强对流个例,地面10 m高风场风速均方根误差降低0.3 m/s,风向均方根误差降低13°;边界层三维风场,风速均方根误差降低0.8 m/s,风向均方根误差降低10°。平原站点融合以后风速、风向预报效果有较大改善,山区站点融合以后改善相对较小。（3）通过对2017年7月20日暴雨和7月7日雷暴大风个例的详细分析,发现融合基于雷达资料四维变分同化获得的高分辨率临近预报风场对各对流系统中的中尺度结构特征给出了更加细致准确的描述。      

线性化物理过程对GRAPES 4DVAR同化的影响

线性化物理过程能够改善四维变分同化中极小化收敛的稳定性和增加极小化过程中对大气物理过程和动力更加精确的描述,它是四维变分同化中非常重要的一部分。通过在GRAPES全球模式中研究线性化物理过程,尤其是两个湿线性化物理过程,改善切线性模式预报精度,来提高GRAPES全球四维变分同化的分析和预报效果。线性化物理过程的开发首先需要简化原非线性化物理过程中的强非线性项,然后对线性化物理过程进行规约化,以抑制切线性扰动的异常增长。目前GRAEPS全球模式中的线性化物理过程主要包括次网格尺度地形参数化、垂直扩散、积云深对流和大尺度凝结。线性化物理过程预报精度的检验方法是通过选择合适大小的初始扰动（同化分析增量）,来比较非线性模式和切线性模式中的扰动演化的纬向平均误差。然后以绝热版本的切线性模式为基础,通过冬、夏两个个例试验来分别检验4个线性化物理过程的12 h预报效果。试验结果表明,通过添加次网格地形参数化和垂直扩散两个干线性化物理过程方案,可以有效抑制住绝热版本切线性模式低层扰动的异常增长,大幅度改善切线性模式预报效果。通过添加积云深对流和大尺度凝结两个湿线性化物理过程,可以在热带区域和中、高纬度地区提高切线性模式中湿变量和温度变量的近似精度,提高GRAPES全球四维变分同化的分析和预报效果。      

NUMERICAL PREDICTION OF AN EXTREME RAINSTORM OVER THE PEARL RIVER DELTA REGION ON 7 MAY 2017 BASED ON WRF-ENKF

An ensemble Kalman filter based on the Weather Research and Forecasting Model (WRF-EnKF) is used to explore the effectiveness of the assimilation of surface observation data in an extreme local rainstorm over the Pearl River Delta region on 7 May 2017. Before the occurrence of rainstorm, the signals of weather forecasts in this case are too weak to be predicted by numerical weather model, but the surface temperature over the urban area are high. The results of this study show that the wind field, temperature, and water vapor are obviously adjusted by assimilating surface data of 10-m wind, 2-m temperature, and 2-m water vapor mixing ratio at 2300 BST 6 May, especially below the height of 2 km. The southerly wind over the Pearl River Delta region is enhanced, and the convergence of wind over the northern Guangzhou city is also enhanced. Additionally, temperature, water vapor mixing ratio and pseudoequivalent potential temperature are obviously increased over the urban region, providing favorable conditions for the occurrence of heavy precipitation. After assimilation, the predictions of 12-h rainfall amount, temperature, and relative humidity are significantly improved, and the rainfall intensity and distribution in this case can be successfully reproduced. Moreover, sensitivity tests suggest that the assimilation of 2-m temperature is the key to predict this extreme rainfall and just assimilating data of surface wind or water vapor is not workable, implying that urban heat island effect may be an important factor in this extreme rainstorm.     

C波段雷达在新疆区域数值预报系统中的应用初探

基于新疆区域数值预报系统（Desert Oasis Gobi Rapid Analysis Forecast System,下称DOGRAFS）,开展了同化C波段雷达资料对2010年10月6日发生在新疆的一次强降水过程预报结果影响的试验分析。其中设计包括不同化任何资料、同化常规资料、同化雷达反射率因子、同化雷达径向风、同时同化反射率因子和径向风五组试验,重点分析了雷达资料同化对此次天气过程降水、温度以及风速模拟效果的影响。结果表明：（1）同化雷达径向风和同时同化径向风和反射率因子相对于其他三组试验,对降水预报的TS和ETs评分更高;（2）相对于其他三组试验,同化雷达径向风和同时同化径向风和反射率因子对模式垂直方向上的温度、风速预报偏差具有一定的改善效果;（3）对于地面2 m温度和10 m风速而言,同化常规观测资料比其他四组试验预报的平均偏差和均方根误差更小;其它四组试验误差相当,差别不明显,表明同化雷达对近地面层温度和风的影响不明显。本研究旨在探索C波段雷达观测资料在新疆区域数值预报系统中的适用性,为今后雷达观测资料在业务系统中的应用提供参考。